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ABSTRACT 

Context. (Ultra)Luminous Infrared Galaxies [(U)LIRGs] are much more numerous at cosmological distances than locally, and are 
likely the precursors of elliptical galaxies. Studies of the physical structure and kinematics of representative samples of these galaxies 
at low redshift are needed in order to understand the interrelated physical processes involved. Integral field spectroscopy (IPS) offers 
the possibility of performing such a detailed analysis. 

Aims. Our goal is to carry out IPS of 42 southern systems which are part of a representative sample of about 70 low redshift (z < 0.26) 
(U)LIRGs, covering different morphologies from spirals to mergers over the entire infrared luminosity range (L//i=10" - IO'^^Lq). 
Methods. The present study is based on optical IPS obtained with the VIMOS instrument on the VLT. 

Results. The first results of the survey are presented here with the study of two galaxies representative of the two major morphological 
types observed in (U)LIRGs, interacting pairs and morphologically- regular, weakly-interacting spirals, respectively. We have found 
that IRAS P06076-2139 consists of two low-intermediate mass (0.15 and 0.4 m.) galaxies with relative velocities of ~ 550 km s"' 
and, therefore, it is unlikely that they will ever merge. The VIMOS IPS has also discovered the presence of a rapidly expanding 
and rotating ring of gas in the southern galaxy (Gj). This ring is interpreted as the result of a nearly central head-on passage of an 
intruder about 140 million years ago. The mass, location and relative velocity of the northern galaxy (G„) rules out this galaxy as 
the intruder IRASP 12115-4656 is a spiral for which we have found a mass of 1.2 m,. The ionized gas shows all the kinematic 
characteristics of a massive (8.7 x 10"'Mo), fast rotating disk. The neutral gas traced by the Nal doublet shows distinct features not 
completely compatible with pure rotation. The neutral and ionized gas components are spatially and kinematically decoupled. The 
analysis presented here illustrates the full potential of IPS in two important aspects: (i) the study of the kinematics and ionization 
structure of complex interacting/merging systems, and (ii) the study of the kinematics of the different gas phases, neutral (cool) and 
ionized (warm), traced by the NaD and Ka lines, respectively. 
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1. Introduction 

Most of the observational efforts devoted in recent years to the 
study of galaxy formation and evolution have been focused on 
obtaining and analyzing large surveys of different characteris- 
tics such as, for instance, 2MASS (Skrutskie et al. 2006), SDSS 
(Adelman-McCarthy et al. 2006), 2DF (Colless et al. 2001), 
UDF (Beckwith et al. 2006), GOODS (Giavalisco et al. 2004), 
COSMOS (Capak et al. 2007), GEMS (Rix et al. 2004), etc. 
Thanks to these studies observational constrains on integrated 
properties such as the galaxy luminosity functions, galaxy (stel- 
lar) masses, sizes, spectral energy distributions, etc., as well as 
their evolution with redshift, have drastically improved (see, for 
instance, Rix et al. 2006 and references therein). However, the 
observational limits imposed by the HST and 10 meter class tele- 
scopes have been reached and a substantial new improvement 
using this methodology will probably have to wait until a new 
generation of telescopes (JWST, ELT, etc) comes into operation. 
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In addition to those integrated properties, a comprehensive 
picture of galaxy formation should explain the internal structure 
of galaxies (i.e., internal kinematics, internal stellar population 
gradients, dust distribution, ionization structure, nuclear prop- 
erties and effects, interaction with the IGM, etc). Such detailed 
studies are more demanding in terms of 'data quality per object' 
requiring high S/N two dimensional spectroscopic information 
with high angular and spectral resolution. Thanks to the advent 
and popularization of Integral Field Spectroscopic instruments 
it is now possible to obtain high quality two-dimensional spec- 
troscopy (i.e. 3D data) over a wide wavelength range (> 1000 A) 
in one spectral setting. Although most such studies have been 
based on local samples of galaxies (e.g., Mediavilla et al. 1997; 
Peletier et al. 1999, 2007; del Burgo et al. 2001; de Zeeuw et 
al. 2002 and references therein), recent work also includes high- 
z studies (e.g.. Smith et al. 2004; Forster-Schreiber et al. 2006; 
Flores et al. 2006; Puech et al. 2007; Law et al. 2007). However, 
limitations in sensitivity and angular resolution make the studies 
of the internal structure of high-z galaxy populations via IFS ex- 
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tremely difficult, and the observation of comprehensive samples 
(i.e., not just individuals at the tip of the luminosity function) 
will require larger aperture telescopes, and higher angular reso- 
lution than those currently available. 

In the context of the high-z universe, the study of representa- 
tive samples of some low-z galaxies populations are of particular 
importance. That is the case of the luminous and ultraluminous 
infrared galaxies (LIRGs, L/R=L[8-1000//ml = 10"-10'2Lo, and 
ULIRGs, L/^>10^^ - IO'^Lq, respectively). These objects are 
thought to be the local counterparts of cosmologically impor- 
tant galaxy populations at z > 1 such as the submillimiter galax- 
ies (SMG; e.g., Frayer et al. 2004), and the majority of the 
infrared sources found in recent Spitzer cosmological surveys 
(e.g., Perez-Gonzalez et al. 2005). ULIRGs and LIRGs have 
large amounts of gas and dust, and are undergoing an intense star 
formation in their (circum)nuclear regions (e.g. Scoville et al. 
2000, Alonso-Herrero et al. 2006, and references therein). This 
starburst activity is believed to be their major energy reservoir, 
although AGN activity may also be present (Genzel et al. 1998). 
In many cases (especially for ULIRGs) these objects show clear 
signs of an on-going merging process producing moderate-mass 
(0.1-1 m,) ellipticals as the end product of the merger (Dasyra 
et al. 2006, Genzel et al. 2001). Although we have a qualitative 
picture, details on how the physical and dynamical processes in- 
volved are interrelated and how they determine the internal struc- 
ture of these cosmologically important objects are far from being 
well understood. 

We have already started a program aimed at studying the 
internal structure and kinematics of a representative sample 
of low-redshift LIRGs and ULIRGs using optical and near-lR 
Integral Field Spectroscopy. These objects are natural laborato- 
ries that allow to investigate in detail the internal structure of 
these galaxy populations in a high S/N, high spatial resolution 
regime, capturing their complex two dimensional structure. We 
are mainly analyzing rest-frame optical and near-lR spectral di- 
agnostic features on hnear scales of about 1 kpc, that will ap- 
pear in the near- and mid-IR in the high-z galaxy populations 
to be studied with future instruments such as the NIRSpec and 
MlRl instruments on board the JWST. This program initially has 
been focused on the most luminous objects, the ULIRGs (e.g., 
CoUna, Arribas & Monreal-Ibero, 2005 and references therein; 
Garcia-Marm 2007) and now is being extended towards lower 
luminosities, i.e. the LIRG luminosity range. 

The present paper is the first of a series of studies based on 
a representative sample of LIRGs observed with VLT-VIMOS. 
With ~ 2000 spectra per object, this study involves a large anal- 
ysis effort. In this first paper we describe the sample, the experi- 
mental set-up used for all the observations with VIMOS, as well 
as the reduction procedures and analysis methods. The analysis 
of two of the galaxies representative of the sample is shown to 
illustrate and emphasize the power of the survey in several areas, 
which after completion will include, among others, the follow- 
ing: (i) two dimensional kinematics of different phases (neutral 
and ionized) of the interstellar gas as a function of luminosity 
and morphology, (ii) two dimensional ionization structure, and 
the interplay between ionization and kinematics (i.e., the role of 
shocks), and (iii) star formation in the (circum)nuclear and ex- 
ternal regions (candidate Tidal Dwarf Galaxies and/or precursors 
of globular clusters) in interacting/merging systems. Throughout 
the paper we will consider Ho=70 km s "'Mpc"' , £Im= 0.7, Qa= 
0.3. 



2. The IFS (U)LIRG Survey and the VIMOS sample 

The IFS (U)LIRG Survey is a large program aimed at the study 
of the two dimensional internal structure of a representative sam- 
ple of LIRGs and ULIRGs using integral field spectroscopic 
facilities both in the southern (VIMOS, Le Fevre et al. 2003), 
and northern hemisphere (INTEGRAL, Arribas et al. 1998 and 
PMAS-CAHA, Roth et al. 2005). The complete survey contains 
about 70 galaxies covering the entire luminous and ultralumi- 
nous infrared luminosity range, and representing the different 
morphologies (spirals, early interactions, advanced mergers, and 
post-mergers) observed in this class of objects. 

The VIMOS sample contains a total of 42 LIRGs and 
ULIRGs (see Table 1 for general properties). The LIRG sam- 
ple consists of 32 systems (42 galaxies) drawn from the IRAS 
Revised Bright Galaxy Sample (RBGS, Sanders et al. 2003), and 
contains all the southern systems from the complete volume- 
limited sample of LIRGs of Alonso-Herrero et al (2006). The 
LIRG subsample has a mean redshift of 0.022 (+ 0.010), with 
distances ranging from ~ 40 to 220 Mpc and covers the entire 
luminosity range. 

The ULIRG sample consists of 10 systems (13 galaxies) se- 
lected from the IRAS 1 Jy sample of ULIRGs (Kim et al. 2002), 
the RBGS (Sanders et al. 2003), and some objects from the 
HSTAVFPC2 snapshot sample (ID 6346 PI: K. Borne), with an 
average redshift of 0.126 (±0.064). 

Galaxies in the VIMOS sample were selected such as to 
cover not only the relevant luminosity range but also represent- 
ing the different morphological types identified in these systems. 
The infrared luminosity, L,r=L(8- lOOOyum), was calculated using 
the fluxes in the four IRAS bands according to the standard pre- 
scription (see Table 1 in Sanders & Mirabel 1996). The values 
obtained are slightly different from those given by Sanders et 
al. (2003) due to the different luminosity distances considered 
(typically our values are ~ 20% higher than those by Sanders et 
al.). Since the definition of LIRGs and ULIRGs is based on a 
luminosity criterion this may induce some apparent inconsisten- 
cies in their classification (for instance, one object could have 
a luminosity higher than lO^^L© [ULIRGs] for us and not for 
Sanders et al.). However, none of the objects in the sample were 
so close to the criterion limit to be in such a situation. 

The morphology class was derived using mainly ground 
based images (Digital Sky Survey available on the NED). For 
several objects these images were not enough to probe the mor- 
phologies, and additional HST images from programs 6346 
(WFPC2) and 10592 (ACS, PI: A. Evans) were used. A sim- 
plified morphological/merging classification scheme similar to 
that given by Veilleux et al. (2002) for ULIRGs has been used. 
Objects that appear to be just single isolated objects without ev- 
idence for a merging process are classified as class 0. Galaxy 
mergers in a pre-coalescence phase, with two well differenti- 
ated nuclei, are classified as class 1, and objects with pecuUar 
morphologies suggesting a post-coalescence merging phase are 
classified as class 2 (see Fig.l and Table 1). This classification is 
based on available optical imaging, and may be revisited as our 
program progresses. In particular, objects classified initially as 
may show characteristics of 2 -like objects after a more detailed 
analysis. Therefore, since at the time of writing some objects in 
the sample are pending observations and reductions, it should be 
considered as indicative only. 

The VIMOS sample is not complete either in flux, lumi- 
nosity or distance. Complete samples of (U)LIGRs, if includ- 
ing the whole luminosity range and morphological types, are so 
numerous that they would require a prohibitive amount of ob- 
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Fig. 1. Distribution of the objects of the sample in the luminosity - redshift (left) and luminosity - morphological/merging type 
(right) planes, where indicates a single isolated object, 1 interacting pairs in a pre-coalecense phase, and 2 a single system with 
evidence of having suffered a merger (i.e. post-coalescence phase). Circles correspond to the current VIMOS sample, whereas 
crosses and filled triangles are the INTEGRAL and PMAS samples, respectively (see text). The dashed Une divides the sample in 
LIRGs (black symbols) and ULIRGs (red symbols). 



serving time. (Moreover, current reduction and analysis tools 
are not able to handle the large volume of data generated by 
very large samples of IPS data.) This is illustrated in Figure 
1, which shows the need to move to higher redshifts to in- 
clude higher luminosity and morphologically diverse objects. As 
mentioned above, in addition to the VIMOS sample, two ad- 
ditional programs in the northern hemisphere using IFS instru- 
ments (INTEGRAL-WHT and PMAS -C AH A) are being carried 
out by us. With INTEGRAL and PMAS we typically observe 
the spectral ranges 4500-7500 A and 3700-7000 A, respectively, 
with a spectral resolution of about 6A in both cases. The spa- 
tial sampling of these data (about 1 arcsec) is somewhat poorer 
than the present VIMOS data. The PMAS sample contains all 
the northern hemisphere targets of the complete volume-limited 
subsample of nearby LIRGs selected from the RBGS (Alonso- 
Herrero et al. 2006). The INTEGRAL sample is a representative 



set of 22 low-z cold and warm ULIRGs covering different mor- 
phological and activity classes (Garcia-Marm 2007). 

In terms of luminosity, the global sample of galaxies 
therefore contains three different subsamples, (i) the complete 
volume-limited local-LIRGs of Alonso-Herrero et al. (2006) (L- 
LIRG sample), (2i) the extended LIRGs (E-LIRG sample), and 
(3i) the ULIRG sample. The distribution of the different sam- 
ples in the L,> - z and L„ - morphology class planes are given in 
Figure 1, where different symbols indicate the instrument used. 

3. Observations, data reduction, and line fitting 

3.1. Observations 

The observations were obtained in service mode during 
semesters P76 and P78 using the integral field unit of VIMOS 
(LeFevre et al. 2003) available at the unit 3 (UT3/Melipal) of 
the Very Large Telescope (VLT). The high resolution mode "HR 
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Orange" (grating GG435), covering approximately the spectral 
range between 5250-7400 A with a spectral resolution of about 
3400 (dispersion of 0.6 A pix"') was used. In this mode the field 
of view covers 27"x27" with a magnification set to 0'.'67 per 
fiber. For this instrument the fibers are closely packed at the 
pseudoslit (mean distance between spectra at the detector is 5 
pixels, with a mean FWHM of 4 pixels). This allows us to use 
40 X 40 fibers obtaining a total of ~ 1600 simultaneous spectra 
per pointing. A square 4 pointing dither pattern with a relative 
off-set of 2'.'7, equivalent to 4 spaxels, was used per galaxy. The 
exposure time was set to 750 (720) s per pointing for a total inte- 
gration of 3000 (2880) seconds per galaxy during cycle 76 (78). 
Weather conditions were in general good, but vary from case to 
case. The mean seeing and air mass for the two systems analysed 
in this paper were 1'.'2, and 1.0 for IRAS F06076-2139 and 1'.'5 
and 1.2 for IRAS F12115-4656, respectively. 

3.2. Data reduction 

The VIMOS data were reduced with a combination of the 
pipeline provided by ESO (version 2.0.2) via "esorex, version 
3.5. r'Q) and a set of costumized IRAF and IDL scripts. Initially 
the pipeline was used to perform the basic reduction steps to the 
individual quadrants (bias, spectra tracing and extraction, cor- 
rection of fiber and pixel transmission, and relative flux cali- 
bration). Once the individual quadrants were partially reduced 
they were combined into a single data-cube. In order to evalu- 
ate the data reduction quality, a series of tests were performed 
after pipeline processing of each individual pointing, and after 
generating the combined data-cube with the four independent 
dithered pointings. Results were similar for the different point- 
ings/objects and are presented in Fig. 2 for one of the pointings 
to illustrate the final quality of the calibrated data. Flux cali- 
bration was checked by comparing the response function ob- 
tained with the ESO pipeline and independently with the IRAF 
tasks standard, sensfunc and calibrate for one of our ob- 
jects. Differences were found to be smaller than 5% between the 
two. The [O i] 6300.3 A sky line was used to check the fiber- 
to-fiber transmission correction and wavelength calibration as 
well as to obtain the instrumental profile (i.e FWHM) for each 
spaxel/spectrum (Fig. 2 illustrates the results obtained). This line 
is the optimal one for these goals since it is bright, isolated, and 
lies between the absorption Nal 5890,5896 A doublet and the 
Hq'-[NII] emission lines complex, which are central to the anal- 
ysis. The effective spectral resolution was established by mea- 
suring the FWHM of a single Gaussian profile fitted to this line 
in each spectrum (~ 1600 in total). The spectral resolution turned 
out to be uniform over the entire field-of-view with an average 
value of 1.83 + 0.25 A (FWHM), equivalent to 87 + 12 km 
s~' at the corresponding wavelength. The obtained mean central 
wavelength (6300.28 + 0.11 A) coincides well with its actual 
value (6300.304 A, Osterbrock et al. 1996) indicating negligible 
systematic errors in velocity (+ 5 km s"' at one sigma level). 
Residuals in flat-fielding affecting the fiber-to-fiber flux calibra- 
tion were derived measuring the relative flux (in counts) of the 
sky line for each of the spaxels in a given sky frame normalized 
to the median flux value of all the spaxels. One-sigma differ- 
ences close to 30% were found, indicating therefore that the rel- 
ative flux of the same emission/absorption line (e.g. He) coming 
from different parts of the same galaxy could have uncertainties 
of up to 30% due to residual calibration effects. The issue of 



flat-fielding and its impact in the relative fiber-to-fiber calibra- 
tion deserves further investigation for the entire set of VIMOS 
data. 

Similar tests to assess the uncertainty in the wavelength cali- 
bration in the blue part of the spectral range were performed with 
the 5577 A sky line. The behaviour of the wavelength calibration 
in the blue end of the spectral range is poor in about half of the 
cases, with errors of up to ~2 A. This is due to the lack of bright 
lines in this spectral range for the available arc calibration lamps 
used during standard operations. For the IFU-HR_orange mode 
the He 50 15. 675 A line is outside the covered spectral range. The 
next blue line (Ne 5400. 56A) is very faint and is not identified 
in most of the cases, being the next usable line at 5862 A. Thus, 
for this VIMOS configuration, the wavelength calibration is well 
constrained only in the ~ 5800-7400 A range. 

Bad spaxels were identified by looking at the output tables 
of the task vmifucalib from the ESO pipeline as well as at the 
extracted continuum lamp spectra. The flat spectra correspond- 
ing to fibers with very low transmission (smaller than 0.25 times 
the median transmission) were rejected. The four pointings were 
combined and the cosmic rays rejected. Because of the dithering 
pattern chosen, most of the information in the dead spaxels was 
recovered during this process. 

Finally the sky was subtracted by obtaining a high S/N ra- 
tio sky spectrum combining a set of spaxels not affected by the 
galaxy emission. In those few cases where the galaxy covered 
the whole VIMOS field of view, the sky spectrum obtained for 
a galaxy observed during the same night was used, scaling its 
flux until the sky line residuals in the subtracted image in the red 
part of the spectrum were negligible. The final "supercube" has 
44x44 spaxels (1936 spectra) . 

3.3. Line Fitting 

The 1936 spectra of each "super-cube" data were inspected vi- 
sually in order to check for possible reduction artifacts, as well 
as for the presence of multiple components, peculiar line pro- 
files, etc. In order to extract the relevant emission line infor- 
mation, line profiles were fitted with a Gaussian model. A sin- 
gle Gaussian profile per emission line was sufficient in most 
of the cases, although occasionally two, or even three compo- 
nents per line were needed. Therefore, in a first step, we fit au- 
tomatically all the lines to a single Gaussian. Then, if necessary, 
multi-component fits were performed in some spaxels restricted 
to certain regions of the field-of-view. These steps were carried 
out with the DIPSO package (Howarth & Murray 1988) and the 
MPFITEXPR algorithm implemented by C.B. Markwardt in the 
IDL environment , which allows us to fix wavelength differ- 
ences and line intensity ratios according to atomic parameters 
when adjusting multiple emission lines such as the Hq'-[Nii] 
spectral range. Several tools have been developed to perform 
these tasks as automatically as possible. By default the same line 
width was considered for all the lines. For each emission line 
we ended with the following information: central wavelength, 
FWHM, and flux intensity. The central wavelengths were trans- 
lated into heliocentric velocities taking into account the radial 
velocities induced by the Earth motion at the time of the obser- 
vation. These were corrected with the IRAF tool rvcorrect. 
The FWHM were corrected spaxel-to-spaxel by the correspond- 
ing instrumental profile (average of 1.83 + 0.25 A, see Fig. 2) 
obtained from the [OI] 6300A sky line. 
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Fig. 2. Histograms showing the results of the calibration checks using the [OI] sky line at 6300.3 A for the first pointing of IRAS 
F06076-2139, after processing the data through the pipeline. Left: line (instrumental) width distribution for all the spaxels, and 
binned in 0.05 A-width intervals. Middle: Distribution of the wavelength for the centroid of the sky Une. Right: Distribution of the 
total line flux normalized to the median flux of the spaxels. 



4. First results and discussion 

In this section we present the first VIMOS results obtained for 
this program. Here we illustrate the full potential of such a sur- 
vey by analyzing two representative galaxies. The whole sample 
will be analysed in a series of future publications. The galaxies 
presented here were selected taking into account their different 
morphological and kinematic characteristics. On the one hand, 
IRAS F06076-2139 is a system with a rather complex morphol- 
ogy showing the presence of two clearly distinct galaxies that 
appear to be interacting. On the other hand, IRAS F12115-4656 
is a spiral galaxy which is in (weak) interaction with a relatively 
distant companion located at about 100 kpc (projected) away. 
Their infrared luminosities are 10^'-^^ and 10^' '^ Lq, respec- 
tively, and therefore they span a factor of 3 .6 in luminosity within 
the range covered by LIRGs. They may be considered represen- 
tative of LIRG sub-types and, as we will show in the following, 
the VIMOS data have revealed interesting new characteristics in 
both cases. The two systems are discussed separately in the rest 
of the paper. 

4.1. IRASF06076-2139 

The ground-based optical (Digital Sky Survey) and near-infrared 
(2MASS) images obtained from the NASA Extragalactic 
Database (NED) for this object show a relatively simple mor- 
phology: a double-nucleus embedded in an elliptical-like galaxy. 
However, the recently available archival ACS/HST images (pro- 
gram 10592, PI: A. Evans) reveal a very rich and complex mor- 
phology of two clearly distinct galaxies (hereafter referred to as 
G„ and G., for the northern and southern galaxy, respectively) 
with a projected nuclear separation of about 9 arcsec (i.e. 6.7 
kpc). The ACS I-band image (i.e. filter F814W) is presented in 
Figure 3 as a reference. The galaxy G„ shows a diffuse double 
ring (ribbon-like), and a spiral of bright distinct knots in the in- 
nermost regions, crossed by a prominent dust lane. On the other 
hand, Gj has a very prominent inner ring of knots of about 0.5" 
(370 pc) in radius (hereafter referred to as the crown, to be dis- 
tinguished from an outer ring - see below) and a more extended 
envelope, elliptical in shape. Of particular interest in the mor- 
phology of galaxy G, is the presence, as inferred from the ACS 
images, of an external ring structure with a semimajor axis of 6 
arcsec (i.e. 4.5 kpc), concentric to the nucleus, and with several 
blue compact knots at different angular locations. Finally, many 
faint blue knots are detected around the main body of the system, 
in particular towards the north and northwest (see Fig. 3). 



4.1 .1 . Spatial structure of the emission line profiles. 
Evidence of the kinematic and ionizing complexity 

The spectra show prominent changes in the Une profiles, velocity 
shifts and emission line ratios, from region to region within the 
system (see Figure 3 for some examples of selected regions). 

The spectra of the two nuclei associated with the two galax- 
ies of the system (N„ and N, in Figure 3, respectively) show 
emission lines with a large relative velocity difference of ~ 550 
km s~\ and with different line profiles. While N„ shows rel- 
atively broad lines (cr= 104 km s"') with a low intensity blue 
wing, the southern nucleus (N^) has narrower lines (cr=63 km 
s~'). There are some other regions (spectra [18,24] and [8,23] in 
Figure 3) where the presence of secondary blue components (i.e. 
blueshifted systems) is clear. Since these regions are spatially 
separated by large distances (about 5 kpc), and show large veloc- 
ity differences, these secondary kinematic components could be 
associated with weak emission from the other galaxy of the sys- 
tem seen in projection along the line of sight. In fact, in certain 
spaxels (spectra [17,26]) two main kinematic systems are iden- 
tified, one per galaxy. Interestingly, the |N1I|/Hq' ratio is quite 
different for these two main velocity systems (see differences 
between the [18,24] and [8,23] spectra). Finally, there are some 
other spectra in the outer regions of the system with narrow lines 
associated with either diffuse gas within the external ring of Gj 
([7,35]), the compact star-forming knot north of G„ (spectrum 
[30,24]), diffuse ionized gas in the outer regions of the system 
(14,5] and [13,7]), or bright blue knots ([10,5]) . Two knots de- 
noted by "RK" outside the VIMOS field of view show very red 
colours (B-I ~ 3.5-4.1), as inferred from the F435W and F814W 
ACS images. 

4.1 .2. Stellar and ionized gas distributions. Detecting new 
structures in the system 

In order to fix the coordinate system for the VIMOS spectral 
maps into the HST reference system, a VIMOS emission line- 
free, stellar continuum image (Figure 4, top-left) was generated 
by integrating the spectra in the observed 6629 - 6733 A spec- 
tral range (rest frame 6390 - 6490 A). Aside from the steps men- 
tioned in Section 3, some spaxels of this image had to be cleaned 
similarly to the standard cosmic ray rejection in standard imag- 
ing. This map follows the global structure of the ACS image 
although, obviously, with poorer spatial resolution. However, 
the patchy structure around the nucleus of the northern galaxy 
and the presence of the crown in the southern nucleus intro- 
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Fig. 3. FULL VERSION at http://damir.iem.csic.es~aalonso/extragalactic/publications/publications.html VIMOS spectra in 
the redshifted Hci'-[NII] spectral range from selected regions of IRAS F06076-2139, as indicated in the HST-ACS F814W-band 
image, which is represented on a logarithmic scale to enhance the low intensity regions. (This image is a subsection of the calibrated 
image as provided by the HST archive, after rotating and removing the cosmic rays using standard IRAF tasks.) The maximum 
(white) and minimum (black) correspond to 6 and 0.1 in arbitrary flux units. Numbers inside the panels identify the spaxel. The 
spectra shown illustrate the changes in velocity, line profiles and emission line ratios. The yellow frame indicates the VIMOS field of 
view represented in other figures of this galaxy in the paper (the actual VIMOS field of view is somewhat larger). The small square 
in the right lower corner of the figure shows the size of a VIMOS spaxel (0.67" x 0.67"). The arrow indicating the orientation has 
a length of 5". The symbol RK identifies the red knots (see text). 



duces some uncertainties in the VIMOS- HST reference sys- 
tems, which have been estimated to be + 0.3". The other panels 
shown in at the top row of Figure 4 were obtained fitting a single 
component per emission line. They correspond to Ha, [Nll]/Ha 
line intensity ratio, mean radial velocity (v), and velocity disper- 
sion (cr). 

The structure of the ionized gas traced by the Ha line looks 
very different from the stellar light distribution mentioned above. 
Galaxy G„ shows a more compact Ha morphology than the stars 
with the Ha emission centered on and dominated by the nucleus 
of the galaxy. On the other hand, Gj exhibits not only emis- 
sion associated with the nucleus (most likely associated with the 
"crown" observed in the ACS image) but also a prominent ring 
of ionized gas of about 16 arcsec in diameter (i.e. 12 kpc) and 
with its major axis along PA ~ 13 degrees. In addition, a chain 
of several Ha clumps forming a bended discontinuous (broken) 
tidal tail-like structure is detected at about 14-16 arcsec (~10.5 
- 12. kpc) west and northwest of the nucleus of G.,. All of these 
Ha clumps, which are not detected in the continuum map, show 
a remarkable positional coincidence with the external blue faint 
knots detected in the ACS F435W image (see Figure 5), indicat- 
ing a clear physical association between them. Therefore, these 
Ha clumps must be young star-forming regions, or tidal dwarf 
galaxy (TDG) candidates, formed as a consequence of a colli- 
sional encounter. 

The observed Ha flux is dominated by the two nuclei of 
the galaxies which contribute 39% (northern nucleus) and 29% 
(southern nucleus) of the total measured flux, respectively. The 
ring accounts for 22.5% of the Ha flux while the rest (9.5%) 
comes from the external chain of Ha clumps. Without a correc- 
tion for internal extinction, the total observed Ha luminosity in 
the system is equivalent to a star formation rate of 1 M© yr"', al- 
most two orders of magnitude lower than the value of 8 1 M© yr"' 
derived from the infrared luminosity (see Table 2 for details). 
Thus, since extinction in (U)LIRGs tends to increase inwards, it 
is likely that the Ha emission is being absorbed preferentially in 
the two nuclei. From the ACS images (see Figures 3 and 5), the 
presence of dust lanes is very prominent in the (circum)nuclear 
regions of the two galaxies. Therefore the observed Ha-derived 
SFR value would reach an agreement with the IR-derived value 
if a visual extinction of about 5 magnitudes exists in at least 
one of the nuclei of the galaxies. Visual extinctions measured in 
the central regions of LIRGs are in the 2 to 6 magnitudes range 
(Alonso-Herrero et al. 2006). 

On the other hand, the total Ha luminosity measured in the 
external Ha clumps corresponds to 1.3 x 10"*" erg s with in- 
dividual clumps (at the VIMOS angular resolution) having Ha 
luminosities in the 0.2 to 0.6 x 10"*° erg s"' range. In addition, 
each of the blue-continuum knots associated with these clumps 
have an effective radius of about 100 to 200 pc (measured di- 
rectly from the ACS B-band image). These properties are simi- 
lar to those measured in some of the TDG candidates detected in 



ULIRGs such as IRAS F08572+3915 and IRAS F15250-1-3609 
at higher redshifts (Monreal-Ibero et al. 2007). 

Finally, the ionization structure of the gas traced by the 
[NII]/l6584/Ha line ratio (Figure 4, top-center panel) is clearly 
distinct between the northern and southern galaxies. The south- 
ern nucleus together with the Ha ring and external clumps all 
have an Hll-like excitation. However, the nucleus and circum- 
nuclear regions of the northern galaxy have ratios changing (in 
log scale) from -0.12 (nucleus) to extreme values of +0.8 at 
about 4 arcsec (3 kpc) from it. Additional line ratios measured 
in this region (log([SII]/Ha) = +0.60) are also high, well above 
the limits of AGN-like galaxies defined by a sample of 85224 
low-redshift emission-line galaxies from the Sloan Digital Sky 
Survey (Kewley et al. 2006). The average spectra in this region 
show a weak Ha in emission on top of a stellar Ha absorp- 
tion line. When corrected for this absorption, the [NII]/Ha and 
[SII]/Ha ratios (in log scale) move to +0.4 and +0.2, respec- 
tively. These ratios correspond to high excitation Seyfert (i.e. 
AGN), LINER-like (weak-AGN or shocks) regions, or a com- 
posite type of ionization. Although a measurement of the [OIII] 
5007A line is needed to distinguish between these two mech- 
anisms, the off-nuclear location of these zones could suggest 
that shocks, as in many other (U)LIRGs (Monreal-Ibero et al. 
2006; Garcfa-Marm 2007), are playing a dominant role in these 
regions. 

4.1.3. Ionized gas kinematics: Towards understanding tine 
nature of the system 

The velocity map (Figure 4, second to the right top panel) shows 
several interesting characteristics that give hints about the nature 
of this system. First, the measured Ha velocity of the nuclei of 
the galaxies (V(Nn)=11186 ± 25 km s"' and V(Ns)=11737 ± 
25 km s"') clearly reveals the presence of two galaxies moving 
with respect to each other with a velocity of about 550 km s 
Since this is a lower limit to their actual relative velocity, it is 
unlikely that these two galaxies form a real merging pair. 

Second, the northern galaxy G„ has an internal velocity gra- 
dient (peak-to-peak) of about 150 km s"' over 6 arcsec (4.5 kpc) 
but with the kinematic major axis tilted by about 60 degrees with 
respect to the photometric axis. The fact that this velocity field 
is consistent with that of a massive, rotating disk of about 5 kpc 
in size and tilted with respect to the stellar mass distribution (see 
§ 4.1.6) would indicate an external origin of the gas. 

Third, the southern galaxy Gj has a regular velocity field 
with the only pecularity that while the kinematic major axis of 
the (circum)nuclear regions (inner 4") is oriented along a PA of 
~ degrees, i.e. close to the stellar major axis of the galaxy, 
the major kinematic axis of the ring is however oriented nearly 
perpendicular to its photometric major axis. This suggests that 
IRAS F06076-2139 is a collisional ring galaxy (e.g. Elmegreen 
& Elmegreen 2006, and references therein). 
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Fig. 4. FULL VERSION at http://dainir.iem.csic.es~aalonso/extragalactic/publications/publications.html VIMOS results for 
IRAS F06076-2139. Upper panels represent the stellar continuum (left) and Ha light (second left) distributions. Changes in the 
ionization conditions are traced by the [NII]/}6584/Ha line ratio (central). The velocity field (second right) and velocity dispersion 
map (right) are also presented. The central and lower rows represent the same characteristics for the two components identified 
in certain regions of the galaxy through a multi-component emission line fit to the observed spectra (see text for details). The 
continuum and Ho- maps are represented on a logarithmic scale between 0.63 (red) and 0.02(black), and 40 (red) and 0.03 (black), 
respectively, in arbitraty flux units. 



Fig. 5. FULL VERSION at http://damir.iem.csic.es~aalonso/extragalactic/publications/publications.htinI Ha emission (con- 
tours) as derived from the VIMOS-IFU data on top of the ACS/HST blue band (filter F435W) image showing the remarkable 
positional coincidence of the external Ha chain of clumps, detected with VIMOS, with the faint, blue stellar blobs located outside 
the main system formed by G„ and Gj. The contours span from 25 to 0.03 in arbitrary flux units in a logarithmic scale. The colors 
range from 10 (white) to 0.03 (black) in arbitrary flux units in a logarithmic scale. 



Finally, the external Ha clumps at distances of about 12 kpc 
from the nucleus of G, shows a rather constant velocity, inde- 
pendent of their orientation and in agreement with the velocities 
measured in the northwestern side of the Ha ring. Thus, this 
large scale velocity structure suggests that the southern galaxy, 
the ring and the external clumps are part of the same galaxy that 
has been involved in an interaction, most likely head-on as indi- 
cated by the velocity structure of the ring (see §4.1.5) 

The velocity dispersion (cr) map (Figure 4, top-right panel) 
also shows a clear distinction between the two galaxies, with G„ 
exhibiting consistently larger dispersions. This suggests that this 
galaxy is more massive than Gj, as also indicated by the light 
distribution. 

Although the northern nucleus (N„) is located in a local max- 
imum of 104 km/s (i.e. it has broader lines than its inner circum- 
nuclear region), it is not well centered on the absolute maximum. 
Actually, the broadest lines are found in the transition region be- 
tween the two galaxies, where two components of similar flux 
but shifted in velocity are present. The velocity dispersion of the 
nucleus of Gj (N^) has a value of 63 km/s, while in the ring and 
chain of clumps is somewhat lower ( 35-40 km/s), indicating a 
more relaxed, quiescent gas. The very large cr values measured 
in a region located ~ 3 arcsec SW from the nucleus of G, are 
due to the the presence of many faint blueshifted components of 
gas, likely connected to the G„. These components are probably 
present along other lines of sight, but they are only noticeable 
when the component of the southern nucleus is faint enough. 

4.1.4. Kinematic deprojection. Origin of tine multiple ionized 
gas velocity components 

As mentioned above, IRAS F06076-2139 shows in some regions 
multiple components in the emission lines. Although a detailed 
analysis of the multiple kinematically distinct components can 
be extremely complex (and may require higher spectral resolu- 
tion), the fact that two main components (in some cases with 
clear evidence of subcomponents) dominate the emission sim- 
plifies this study. Therefore, in the cases where only one sin- 
gle Gaussian per emission Une was not a good representation 
of the observed line profile, two Gaussians per line were fitted 
instead. Note that the spatial identification of the components 
(kinematic deprojection, see Arribas et al. 1996) can be done 
unambiguously thanks to the fact that they have quite different 
kinematic properties (i.e. mean radial velocity and velocity dis- 
persion). Therefore, in general, at each spatial position there was 
no doubt as to which fit should be associated with which compo- 
nent. This allows us to construct spectral maps of each kinemati- 
cally distinct component. The central and lower rows in Figure 4 



represent the spatial distribution of Ha, the [NII]/Ha line inten- 
sity ratio, mean radial velocity (v), and velocity dispersion (cr) 
for the two components. The double components are found only 
in two regions: i) the (projected) interface region between both 
galaxies, and ii) westwards of the southern nucleus. A third re- 
gion with double components was also found to the north of G„ 
but it has not been considered here due to the large uncertainties 
in the fits with two components. 

Looking at these spatial distributions the so called compo- 
nents 1 and 2 share the kinematic and ionization properties of 
galaxies G„ and Gj, respectively, and are therefore physically 
associated with each of them. In fact component 2 (Figure 4, 
bottom panels) has the line widths, mean radial velocities, and 
[NII]/Ha ratios expected for G, at the spatial location where it is 
found. Similarly, the properties of component 1 (Figure 4, cen- 
tral panels) are those expected for G„. However, interestingly, 
the radial velocities of the region close to the southern nucleus 
are intermediate between those expected for the two galaxies. 

4.1 .5. Kinematic evidence for an expanding collisional ring in 

One of the most interesting kinematic features of the velocity 
field associated with G, is the regular pattern of the radial ve- 
locities in the Ha ring. The fact that the largest velocity gradient 
is nearly perpendicular to the photometric major axis of the ring 
suggests that the ring may have a peculiar velocity pattern with 
an important expansion component. 

As a first step we selected the spaxels/spectra containing in- 
formation about the ring. For those spectra showing two compo- 
nents per emission line, only the one corresponding to the ring 
was considered (see previous subsection). To analyze in detail 
the velocity field of the ring, a simple kinematic model (Mihalas 
& Binney, 1981) was fitted to the data: 

v,.(R, 9)^ [n(/?, 9) X sini9-9„) + Q(/?, 9) x cosi9-9„)] x sini + 
Z{R,9)xcosi + Vsys,i^) 

where n(R, 9), Q.{R, 9), and Z{R, 9) are radial, azimuthal, and 
perpendicular components of the velocity field, / is the inclina- 
tion angle (defined as the angle between the normal to the galaxy 
disk [ring] and the line of sight), Vsxs is the systemic velocity, 
9,, the position angle of the projected major axis of the galaxy 
(ring), R the galactocentric distance , and 9 the galactic azimuthal 
angle. This model is rather generic, but it can be simplified for 
the present case taking into account some observational facts. 
First, our data indicate azimuthal symmetry and therefore the 
amplitude of the of the velocity components does not depend 
on the galactic azimuthal angle. Second, the global systemic ve- 
locity of the ring is in agreement with the radial velocity of the 
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Fig. 6. Top. Best fit (red line) and data points for the kinematic 
model of the ring in G, as described in equation (2). Bottom. 
Decomposition of the total radial velocity of the modeled ring 
(red line) into its orbital (green) and azimuthal (blue) compo- 
nents. The dashed line shows the systemic velocity of the ring. 

nuclear region (see Table 2), indicating that the value of the per- 
pendicular component (i.e. Z) is negligible. Taking into account 
these facts, equation (1) can be reduced to: 

VriO) = [n{R)sinie - 0o) + Q.{R)cos{e - Bo)\ x sini + V,,,,, (2) 
where Qo can be easily inferred assuming that the actual geom- 
etry is circular. The observed position angle inferred from these 
images is 1 3" + 2. The measured values for the lengths of the ma- 
jor and minor axes of the observed ring are then used to derive 
the angle of inclination, /, which is found to be 57" + 2 (with a 
global uncertainty of 180"). 

Considering these values, equation (2) is solved for IT, Q, 
and Vjvs by applying a least square fit to the data (see Figure 6, 
top panel, for the 206 data points and the best fit). The best fit 
(Figure 6, bottom panel) was found for an expanding velocity H 
= 70 + 3 kms^\ a clockwise rotation Q = 61+ 3 kms^^, and a 
systemic velocity Vsys = 1 1750 + 1 kms^^ . For H and Q the er- 
rors were estimated from the uncertainties associated with the 9o, 
and / only. For Vsys it was estimated from the standard deviation 
of the differences data - model, which was 0.7 kms^^ . All these 
errors do not include systematic effects due to both the measured 
data, and the model approach. Mathematically the radial com- 
ponent, n, could be explained by both expansion or contraction 
(due to 180" uncertainty in /). However, from a physical point of 
view, we can safely exclude radial contraction at such a velocity. 
Furthermore, in ring galaxies we actually expect the presence of 
expansion induced by the wave created in reaction to the impact. 
The fact that the ring is expanding fixes its own geometry in 
such a way that its SW sector must be the one closer to us. That 
in turn fixes the direction of the azimuthal (rotational) compo- 
nent which is clockwise. This result, which relies on solid kine- 
matical arguments, seems to contradict the morphological spiral 
appearance of some structures inside the ring in the ACS B-band 



image, which suggest an anti-clockwise rotation if it is assumed 
that these structures are trailing. However, this apparent contra- 
diction can be explained as differential rotation inside the ring 
producing a rising 'rotation curve' such that the outer parts are 
rotating faster than expected for solid-body rotation. 

4.1 .6. Evidence for a gaseous kinematically decoupled 
rotating disk in G„ 

The velocity field of G„ does not show a clear large-scale rota- 
tional pattern. The kinematic major axis seems to be tilted with 
respect to the stellar photometric axis, and the velocity disper- 
sion map is also difficult to interpret as due to pure rotation. The 
ACS images also shows a quite complex stellar light distribution 
with several dust lanes and what appears to be star-forming re- 
gions in a spiral-like structure (Figure 3). We attempt to fit the 
velocity field to the model described above (eq.[3]). As a first 
order approximation, we considered a rotational model with a 
constant value for Q. The best fit model for /=58 degrees (de- 
rived from the ACS images) and considering Q, 0o, and Vsys as 
free parameters gave a rotational velocity of 76 km s"', a po- 
sition angle for the major axis of +40 degrees, and a systemic 
velocity of 11236 with an rms of ~ 3 km s"'. If Og is fixed to 
-20 degrees (i.e. the stellar photometric axis position as derived 
from the ACS images) the fit has residuals substantially larger 
than the model above. These results indicate that, despite the 
relatively large deviations, the velocity field in galaxy G„ is con- 
sistent with rotation around an axis tilted by about 60 degrees 
from the stellar major axis. This suggests that the ionized gas 
is decoupled from the stars. Moreover, the high velocity disper- 
sions outside the nuclear regions also indicate that the ionized 
gas in G„ has quite complex kinematics. The distribution and 
kinematics of the ionized gas and the stars are often known to be 
decoupled in early-type galaxies (e.g., Bertola & Corsini 1999). 
The velocity fields show misalignments between the stars and 
the gas (e.g., Kannappan & Fabricant 2001; Sarzi et al. 2006). 
This has been used to constrain the origin of the gas in early 
type galaxies. In the case of IRAS F06076-2139, which clearly 
shows strong merging/interaction features, it is not surprising to 
find such misaligment. 

4.1.7. IRAS F06076-2139. Overall kinematic model and 
physical parameters 

The expanding rotational ring model for the southern galaxy has 
been combined with the tilted rotational disk for the northern 
galaxy in order to compare the overall model with the observed 
velocity field of the system (see Figure 7). The model does not 
include the internal regions of the southern galaxy due to the 
limited angular resolution and to the fact that the ACS images 
show the presence of a nuclear star-forming ring (the crown) that 
could have its own velocity field independent of the one detected 
in the outer expanding ring. The model explains well the overall 
velocity field on large scales as the average residual (observed - 
model) is -4 ± 1 1 km s"' and +0.6 ± 29 km s"' for G., (ring) and 
G„, respectively. However, significant residuals, blueshifts larger 
than - 50 km s ' with respect to rotation, appear on smaller 
scales in G„. The nature of these shifts is unclear, but certainly 
demonstrate that the velocity field in G„ is not only due to the 
pure rotation of a tilted disk but is also affected by additional 
radial flows, or tidally-induced motions. 

From the azimuthal velocity component (i.e. rotation), O, we 
can estimate the mass inner to the ring (i.e. 10.1 kpc) in G,, at a 
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Fig. 7. FULL VERSION at http://damir.iem.csic.es~aalonso/extragalactic/publications/publications.html Comparison of the 
Ha observed (left) and modeled (centre) velocity field for IRAS F06076-2139. The modelled velocity field combines the model of 
the tilted rotating disk (case 1, see text) in G„ with the expanding, rotating ring in G,. The residual (observed - modeled) velocity 
field is also given (right). The central spaxel at (0,0) corresponds to the nucleus of G„, and has been excluded since it is a singularity 
for the model. 



radius of 5.6 kpc for G„. The dynamical masses inside these radii 
are M(G,)= 1.7 x 1O"'M0 and M(G„)= 7.5 x 1O'*M0, respec- 
tively. However, these masses represent however only a fraction 
of the overall mass of the system. Under virialization and con- 
sidering a polytropic mass distribution, the overall mass of each 
of the galaxies is given as 1.24 x 10"^ x R/„„(kpc) x cr^(km s"') 
(see Colina et al. 2005 for details) in units of an m, galaxy (1.4 
X 10"MO, Cole et al. 2001). Assuming that Hght traces mass, 
the half-mass radius can be obtained measuring the effective ra- 
dius. For this purpose, polygonal apertures equivalent to circular 
apertures of 12.8" and 8.6" in radius for G„ and G, respectively, 
have been used in the recently available ACS/F8 14W image. The 
polygonal apertures allow us to better estimate the effective ra- 
dius of each of the galaxies as the aperture is selected in such 
a way as to minimize the light contribution from the compan- 
ion. The corresponding effective radii are 2.76 kpc and 2.73 kpc 
for G„ and G,, respectively. Taking into account the measured 
central velocity dispersions of 104 km s"' and 63 km s IRAS 
F06076-2139 is formed by two low-intermediate mass galaxies 
with masses of 5.2 (G„) and 1.9 (G,) x lO"' Mq, i.e. 0.4 and 0.15 
m*, respectively. 

Finally, the expansion velocity of the ring can be used to date 
the epoch when the impact that generated the ring took place. 
Assuming that the ring has been expanding at a constant speed 
of 70 km s"' and that the distance to the nucleus is 10. 1 kpc, such 
an impact happened about 1 .4 x 10*^ years ago. By measuring the 
"width" of the ring, one can also estimate the spread in age of its 
star forming regions, which is ~ 2 x 10^ years. 

4.1.8. Status and evolution of IRAS F06076-2139 

Based purely on the morphology as seen in the high angular res- 
olution ACS images, IRAS F06076-2139 would be considered 
a merging system formed by two bright spirals. However, the 
VIMOS IFS provides a set of new quantitative data that chal- 
lenges this simple scenario, and shows the full complexity of the 
system. 

On the one hand, the newly identified Ho- ring in Gj has been 
interpreted as an expanding rotating collisional ring. Such rings 
are generally produced by head-on collisions with an intruder 
having a mass lower than the victim. Additional morphological 
and kinematic features of the ring give new information about 
this collision. The collision happened about 140 million years 
ago as derived by the size and expansion velocity of the ring. 
The fact that the nucleus of G., lies close to the center of the 
ring indicates that the passage of the intruder was near-central. 
Moreover, the expansion velocity of the ring is relatively high 
(70 km s" ' ) and about 1.15 times its rotational velocity. Since the 
amplitude of the expansion velocity is directly proportional to 
the mass of the intruder (Struck-Marcell & Appleton 1987), the 
intruder should have had a mass smaller than but close to that of 
G^. These new facts rule out a direct collision between G„ and G, 
as the cause of the ring since: i) the intruder (remnants) should 
be located along the minor axis of the ring (G« is more than 
70 degrees off from that axis), and ii) G„ is considerably more 
massive than Gj (by a factor of 2-3), according to the different 
tracers used. The brightest of the blue external knots identified 



in the ACS B-band image, and also detected as an Ha emitter 
with VIMOS, is located at about the right distance and along the 
minor axis of the expanding ring, as expected for the intruder 
(see Figure 5). However, this knot appears to be a member of 
the chain of several clumps detected north and northwest of G., 
sharing similar kinematics and ionization properties and without 
evidence of an old stellar population (traced by the I-band light) 
that could be considered the renmant of the intruder. However, 
the red knots found by comparing the F435W and F814W ACS 
images (see Fig. 3, labels RK), with B-I colors of 3.5 and 4. 1, are 
good candidates for the remnants of the intruder. Unfortunately 
they are outside the VIMOS field of view, and therefore we lack 
the kinematic information to check this hypothesis. 

As for the fate of this system, the high relative velocity be- 
tween the two galaxies is larger than the escape velocity consid- 
ering the galaxy masses infeiTed from the present observations. 
Therefore, it is very unlikely that these two galaxies will ever 
merge, although they could pass through each another. Without 
knowing their actual distance, nor their relative velocity in the 
plane of the sky it is difficult to predict the future evolution with 
certainty. 

4.2. IRAS F1 211 5-4656 

This galaxy is member of a pair included in the catalogue 
of interacting galaxies of Arp & Madore (1987) with the 
second member of the system (IRAS F121 12-4659) located 
at a distance of about four arcmin (89 kpc) SW of IRAS 
F121 15-4656. The emission line spectrum of the nucleus has 
been classified as Hll-like, and its absorption spectrum as due 
to a young stellar population (Sekiguchi & Wolstencroft 1992). 
IRAS F121 15-4656 was also included in the study of Palunas 
and Williams (2000), who presented axisymmetric maximum 
disk mass models for a large sample of spiral galaxies based 
on rotation curves extracted from their two-dimensional Fabry- 
Perot Ha data (Schommer et al. 1993). For the present paper, 
which presents the first results from the VIMOS sample, this 
galaxy was selected mainly because it is representative of a class 
of local LIRGs with no clear morphological indication of being 
involved in a large gravitational interaction or merging process. 
Unlike IRAS F0607-2139, IRAS F121 15-4656 is classified as 
an early type spiral (SA(rs)ab, Buta 1996). 

4.2.1 . Stellar and ionized gas distributions. Tracing the 
ionization sources 

The stellar light distribution was obtained generating a narrow 
emission line-free continuum image in the observed 6508 - 6610 
A spectral range (Figure 8, left panel). The bright knot at about 
10 arcsec NE from the nucleus is a field star. The emission lines 
in the spectra of IRAS F121 15-4656 show in general simple 
profiles with only one kinematic component per line. However, 
some spectra have a faint red (blue) wing predominantly in the 
NW (SB) of the nucleus. In a few cases we found some hints of 
double components, but the relative velocity between these com- 
ponents is at the limit of being resolved with the resolution and 
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S/N level of the present data. Therefore, all the spectra were sat- 
isfactorily fitted with a single Gaussian per emission Une. From 
these fits the Hor light distribution is obtained (see Figure 8, sec- 
ond left panel). 

The general morphology of the ionized gas (traced by the Ha 
map) follows that of the stars. However there are some clear dif- 
ferences. The stellar distribution shows an elHptical shape with a 
well defined central nucleus which appears as a relatively weak 
Ha emitter. The circumnuclear regions show, on the other hand, 
the presence of compact, bright Ha knots. These knots at dis- 
tances of 1 to 1 .5 kpc are most likely HII regions located in the 
irmer spiral arms of the galaxy. Also two regions at about 3.7 kpc 
southwest of the nucleus, and undetected in the continuum, are 
clearly distinct in Ha. Finally, two faint Ha emitting regions at 
5.1 kpc NW and SE the nucleus are also visible. 

The location and nature of the ionizing sources are traced by 
the [NII]/Ha line ratio (Figure 8, central panel). According to 
their ionization status, three main, distinct regions can be iden- 
tified: (a) the Ha-weak nucleus, (b) the circumnuclear Hll re- 
gions, and (c) the elliptical outer ring of enhanced [Nil] emis- 
sion. The weak-Hor nucleus has a relatively high [NII]/Hq; value 
(log[NII]/Ha ~ -0.15). This ratio together with other ioniza- 
tion tracers (log[OI]/Ha of -1.31 and log[SlI]/Ha ~ -0.39) lo- 
cates the nucleus in the borderhne region of Seyfert/LINERs and 
HII regions, according to the standard emission Une diagrams 
(Kewley et al. 2006). Additional [OIII]/I^ line ratios would be 
required to establish a more precise classification. The circum- 
nuclear Ha knots, and the faint two regions located SW of the 
nucleus have line ratios (log[NII]/Ha of -0.37, log[OI]/Ha of 
-1.67 and log[SII]/Ha ~ -0.55) typical of HE regions. Finally, 
the elliptical outer ring clearly defines the boundary of the main 
body of the galaxy in the Ha map. In this elliptical ring the aver- 
age log[Nll]/Ha ratio is -0.02, typical of LINER or Seyfert-like 
ionization. The same elUptical structure is also identified in the 
velocity dispersion map (Fig. 8 right panel) indicating that there 
is a clear spatial correlation between the ionization status of the 
gas and its kinematics. Similar correlations have been found in 
the outer parts of ULIRGs, being explained by the presence of 
local ionizing sources (shocks) in these regions, rather than by 
the ionization due to a central, weak AGN (Monreal-Ibero et al. 
2006). 

Assuming that the gas is ionized only by young stars, the 
equivalent star formation rate associated with the observed Ha 
flux corresponds to 1.9 Mq yr"\ or about 10% of the total 
SFR derived from the IR luminosity (see Table 2 for values). 
Therefore, in order to bring the two values to an agreement, ex- 
tinctions of about 3-4 visual magnitudes in the nuclear regions 
are required. These values are within the range of extinctions 
typical of LIRGs (Alonso-Herrero et al. 2006). 

4.2.2. Ionized gas kinematics. IVIodeling the velocity field and 
measuring tine dynamical mass 

The velocity field of the ionized gas exhibits a regular veloc- 
ity field, typical of rotational motions. The major kinematic axis 
is well aligned with the stellar and ionized gas distributions 
traced by the optical continuum and Ha maps, respectively (see 
Figure 8). In addition, the velocity dispersion has its maximum 
value at the stellar center of the galaxy (within the angular res- 
olution), as expected in galaxies dominated by pure rotation. 
Therefore, the observed velocity field has been modeled with 
the kinematic model described in equation (1) assuming that 
the radial and perpendicular velocity components are negligible 
(YliR, 0) = Z(R, ff) = 0). Considering an inclination value of 55 
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Fig. 9. Rotation curves for the ionized (soUd Une) and neutral 
gas (dashed line). These curves represent the variation of the az- 
imuthal component of the velocity field, Q, with the galactocen- 
tric distance, as inferred after fitting the data to a rotating disk 
model (see text). The nominal uncertainty (rms) in radial veloc- 
ity at each of the evaluated radii is typically about 1 km s^', 
except for the inner ( < 0.4 kpc) and outer (> 6.5 kpc) regions, 
which is somewhat larger (3-5 km s~^). 



degrees (inferred from the stellar continuum map) the 908 data 
points were divided into 12 rings of spaxels located at increas- 
ing galactocentric distances and then equation (1) was solved 
with n, Vsys, and 0o as free parameters. The standard deviation 
was typicaUy lower than 1 km s~' for most of the 12 indepen- 
dent fits performed, and only reached values of a few km s~^ for 
the inner ( < 0.4 kpc) and outer (> 6.5 kpc) regions. Given that 
the values at large galactocentric distances were obtained using a 
limited range in azimuth they should be less reliable and are not 
considered further here. The variation of the azimuthal velocity 
component, Q(R), is presented in Figure 9. This rotation curve 
reaches its maximun value (261 km s"') at a galactocentric dis- 
tance of ~ 3.5 kpc, then it keeps a relatively constant value of 
254 + 4 km s"' up to 6 kpc, when it starts to decline. The mean 
values for the systemic velocity (Vsys) and the PA of the major 
kinematic axis (6o) evaluated at r < 6.5 kpc were 5449 + 11 km 
s"', and 115 + 1 degrees (with an additional uncertainty of ±180 
degrees due to the uncertainty in i), respectively. The results of 
the fit agree very well with previous values (Q= 262km s~\ and 
00 = 294) obtained by Palunas & Williams (2000) from Fabry- 
Perot Ha maps. However, the heliocentric velocity derived from 
the kinematic model disagrees with that given by Strauss et al 
(1992) who found 5543 ± 35 km s~\ but shows a good coin- 
cidence with the value provided by Sekiguchi et al. (1992) of 
5436 ± 38 km s"^ If the spiral arms seen in optical images of 
this galaxy (e.g. in the DSS images) are trailing, the observed 
kinematics indicates anti-clockwise rotation, and the SW sector 
is the closest to us. 

The observed and best-fit modeled velocity field together 
with the corresponding two-dimensional map velocity residuals 
(observed minus modeled) are presented in Figure 10. In gen- 
eral, the rotating disk model is a very good representation of the 
main motions of IRAS F12115-4656, with an average residual 
of 0.2 +12 km s~^ (note that the contiguous E and W spaxels to 
the nucleus have relatively large residuals of about 40 km s"' to 
the red and blue, respectively, due to the Umited spatial resolu- 
tion of the data). However, significant low ampUtude (~ 20 -30 
km s"') velocity residuals along PA 25 degrees are identified 
(see Figure 10, right-upper panel and §4.2.4 for a more detailed 
discussion). 

The dynamical mass of the galaxy within a radius of 6.0 kpc 
corresponds to 8.7 x lO'" Mq for a derived rotation velocity of 

251+2 km s"' (see Figure 9). Note however, that if the velocity 
field of the neutral gas is considered (see next section), the dy- 
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Fig. 8. FULL VERSION at http://dainir.iem.csic.es~aalonso/extragalactic/publications/publications.html VIMOS results for 
IRAS F121 15-4656 showing the stellar continuum (left), and the warm ionized (second left) light distributions together with the 
ionization map traced by the ratio of the [Nil] to the Ha lines (center). Also shown are the velocity field (second right) and velocity 
dispersion map (right) as obtained from the fit of the Ha line to one component Gaussian profile. The continuum and Ha maps are 
represented in a logarithmic scale between 12.5 (red) and 0.16 (black), and 40 (red) and 0.03 (black), respectively, in arbitraty flux 
units. 



Fig. 10. FULL VERSION at http://damir.iem.csic.es~aalonso/extragalactic/publications/publications.htmI Figure presenting 
the results for the velocity fields of the ionized (top panels) and neutral (bottom panels) gas in IRAS F121 15-4656. Observed (left 
panel), modelled (central panel) and residuals (observed minus modelled; right panel) are shown (see text for details). The map 
corresponding to the neutral gas, which was inferred from the Nal lines, has been smoothed by averaging the spectra around each 
spaxel in order to increase the S/N. 



namical mass within the same radius will be about a factor of 2 
lower. As for IRAS F06076-2139 (see §4.1.7), a measurement 
of the total mass of the galaxy can be derived from the value of 
the central velocity dispersion and the effective radius. Since the 
VIMOS field-of-view is so small compared to the diameter of the 
galaxy, the VIMOS continuum image cannot be used to establish 
the effective radius of the light distribution. Instead, the effective 
radius measured in the ESO-Uppsala Galaxy Survey (Lauberts 
& Valentijn 1989) in the B-band light is used. For an eflFective 
radius of 20.6 arcsec (7.7 kpc) and a central velocity dispersion 
of 111 km s \ a total mass of 1.65 X 10'^ Mq, equivalent to an 
1.2 m, galaxy, is derived. 

4.2.3. Interstellar gas phases. Evidence for a kinematic and 
spatial decoupling of the neutral and ionized gas 
components 

The physical properties of the cool neutral interstellar gas in 
galaxies can be obtained by measuring the characteristics of in- 
terstellar absorption lines such as the Na I doublet 5890, 5896 A 
in the optical (e.g. Martin, 2006; Vaisanen et al. 2007). However, 
little is known so far about the two-dimensional kinematics of 
the neutral gas phase in (U)LIRGs due to the difficulty of detect- 
ing these lines. Unlike the emission lines with large equivalent 
widths of typically 10 to 100 A, the Nal Unes have small equiv- 
alent widths of a few A . The high quality of the VIMOS data 
now allows us to measure the spatial distribution and kinematics 
of the neutral interstellar gas in two dimensions, making possible 
a direct one-to-one comparison with those derived for the warm 
ionized gas (see Figure 11). 

There are a number of remarkable features both in the spa- 
tial distribution and in the kinematics of the neutral gas that re- 
quire further discussion. The neutral gas distribution, traced by 
the equivalent width of the Na I 5890 A line, is similar to that of 
the stars with a clear peak at the photometric center of the galaxy, 
and extended emission along the same major axis as the stellar 
distribution (Figure 11). The ionized gas, on the other hand, al- 
though sharing the same major axis, shows a very different spa- 
tial distribution dominated by the circumnuclear regions where 
the young ionizing star forming regions are located (see §4.2.1). 
Moreover, the fact that the equivalent width of Ha in the center 
of the galaxy is low, exactly where Na I has its maximum, could 
partly be due to an increase in the extinction towards the nuclear 
region where a larger concentration of cold dust associated with 
the cool neutral gas would exist. 

When compared with the kinematics of the ionized gas, both 
the velocity field and the velocity dispersion map of the neu- 
tral gas show distinct features (see Figure 1 1). On the one hand. 



the neutral gas shows indications of being affected by rotation 
as traced by the large peak-to-peak velocity amplitude (380 km 
s"') along a major axis centered in the nucleus of the galaxy. 
However, the velocity field shows some departures from a regu- 
lar rotation pattern, and the maximum of the velocity dispersion 
is not at the dynamical center as expected in typical rotating sys- 
tems. Moreover, the velocity dispersion of the neutral gas shows 
very high values (cr ~ 180 km"') at distances of about 0.8 kpc 
north and 2.7 kpc northwest of the nucleus, much higher than 
those measured in the ionized gas (cr ~ 30 to 40 km"') at the 
same locations and with no clear correspondence with any of 
the circumnuclear star-forming regions. Broad Na I lines such as 
those identified in the northwest regions of IRAS F121 15-4656 
have also been detected in some ULIRGs outside the nucleus 
(Martin 2006). 

To further compare the neutral and ionized gas kinematics, 
the rotational velocity pattern of the neutral gas has been sub- 
stracted from that of the ionized gas. This map (Figure 11, cen- 
tral lower panel) shows the morphology of a bipolar flow with 
velocities of up to + 100 km s"', and with its axis oriented 
close to the major kinematic axis of the ionized gas. This be- 
haviour can be explained by the neutral gas rotating slower than 
the ionized gas but along a similar major kinematic axis. This 
scenario is confirmed when modeling the neutral gas velocity 
field with a rotating disk structure in a similar manner as for 
the ionized gas (see Figure 10 lower panels) . As can be seen 
in these panels the model is also a good representation of the 
observed velocity field, with average residuals of 0.1 + 17 km 
s"'. The average systemic velocity at distances of less than 5.5 
kpc is 5460 ±18 km s"' , in agreement with that of the ionized 
gas within the uncertainties. However, the rotation curve of the 
neutral gas (see Figure 9) indicates a rotational velocity system- 
atically lower than that of the ionized gas by ~ 70 - 90 km s"'. 
Also, the kinematic major axis is oriented along PA 98 + 3 de- 
grees, significantly misaligned with respect to that of the ionized 
gas by 17 + 3 degrees. 

Regarding the velocity dispersions, the central values for 
both neutral (116 km s"' ) and ionized (111 km s"') gas agree 
well, as expected if these values trace the mass in the nuclear 
regions. However, outside the nucleus, the overall residuals in- 
dicate that the velocity dispersion of the neutral gas is always 
higher than that of the ionized gas by 35 to 150 km s"' (see 
Figure 1 1 , lower right panel for a two-dimensional map of the 
differences). 

The V/cr ratio for the neutral and ionized gas shows even 
clearer differences. While the V/cr values of the ionized gas are 
high (typically > 3) indicating a system dominated by rotation as 
in spirals, the corresponding V/cr values for the neutral gas are 



12 



(VIMOS (U)LIRGs first results: VLT-VIMOS integral field spectroscopy of Luminous and Ultraluminous Infrared Galaxies 



Fig. 11. FULL VERSION at http://dainir.iem.csic.es~aalonso/extragalactic/publications/publications.html Spatial distribu- 
tion and kinematic properties of the neutral gas in IRAS F121 15-4656 as traced by the Nal doublet (top panels). The equivalent 
width of the Ho- line tracing the regions of ionized gas associated with young star-forming regions is given for comparison (bottom 
left panel). The residuals of the velocity field and velocity dispersion of the neutral gas after subtraction of those measured for the 
ionized gas are also given (central and right bottom panels) to highlight the differences in the kinematic properties of the two gas 
phases. 



lower (typically < 1 ) suggesting the presence of a hotter system 
more affected by the random motions, like the stars in ellipticals. 

Summarizing, the kinematics of the neutral and ionized gas 
demonstrates that the dynamical status of the two gas phases is 
different, one (ionized gas) traces the presence of a rotating star- 
forming disk in the galaxy while the other (neutral gas) could be 
tracing the overall dynamics of the stellar distribution with some 
departures in specific regions. 

4.2.4. Evidence for non-rotational motion of tine ionized and 
neutral gases: correlation with the velocity dispersion 

As already mentioned (§4.2.2) significant low amplitude (~ 20 
-30 km s"' ) velocity residuals along ~ PA -25 degrees have been 
identified (see Figure 10, right-upper panel). Considering the 
kinematic model given by equation (1), such departures from 
pure rotation can be interpreted as radial motion in the plane 
of the disc (i.e., n{R,6) ^ in eq. [1] ) and/or as motion per- 
pendicular to the plane of the disc (i.e., Z(R,0) in eq. [1] ). 
Purely radial motion produces a pattern such that regions sym- 
metric with respect to the nucleus have opposite signs in their 
velocities. This is not observed, thus favoring the hypothesis 
that they mainly represent motions perpendicular to the plane. 
Furthermore, taking into account the geometry of the disc, with 
the SW sector closest to us, the residuals suggest motions out- 
wards the disc. Interestingly, these regions with non-rotational 
motions do spatially correlate with regions of increased veloc- 
ity dispersion (see right panel in Figure 8), indicating a physical 
connection between the two kinematic tracers. 

The same kind of behaviour, and even clearly detected, is 
found in the neutral gas. The regions where the largest departures 
from pure rotation are the detected (up to 50 km s"') also exhibit 
the largest velocity dispersions. In general the velocity residual 
map (Figure 10 lower-right) has a remarkable spatial correlation 
with the velocity dispersion derived from the Nal lines (Figure 
11 upper-right). A correlation between the measured velocity 
dispersions and the presence of non-rotational motion (mainly 
produced by tidal forces) has been reported at a higher inter- 
action/merger scale in ULlRGs (Colina et al. 2005). Although 
IRAS F121 15-4656 has only a weak interaction, we could be 
witnessing a low-scale effect of the same phenomenon. 

5. Conclusions 

This paper presents the first results obtained from the inte- 
gral field spectroscopic survey of local luminous and ultralu- 
minous infrared galaxies being carried out with the VIMOS in- 
strument on the VLT. This sample includes 42 local (U)LIRGs 
systems (65 galaxies) and it is part of a larger representative 
sample of about 70 systems (z < 0.26) also observed from the 
northern hemisphere, which covers the different morphologies 
from spirals to mergers and spans the entire luminosity range 
Lm=10" - IQ^^-^Lq. The two galaxies analysed here (IRAS 
F06076-2139 and IRAS F121 15-4656) have been selected as 
prototypes of the most common morphological types among 



(U)LIRGs, i.e. interacting pairs and regular/weakly-interacting 
spirals, while also covering a wide range in IR luminosity. 

A) The VIMOS data of IRAS F06076-2139 clearly illus- 
trate the potential of IFS in disentangling the complex ionizing 
and kinematic structure detected in pairs or mergers. The main 
conclusions for this system are: 

al) We have found that this system is formed by two inde- 
pendent low-intermediate mass galaxies, the G„ (northern) and 
Gj(southern), of 0.4 m» and 0.15 m», respectively, and with a 
relative projected velocity of about 550 km s"'. Therefore, it is 
unlikely that these galaxies will ever merge. G„ is characterized 
by an AGN-like ionization, and large velocity dispersions (> 100 
km s '), while G, has HII characteristics and lower velocity dis- 
persions (by at least a factor of two) than G„. 

a2) The morphology and velocity field of the ionized gas in 
Gs is consistent with the presence of an expanding, rotating ring 
of about 20 kpc (deprojected major axis) in size. This expand- 
ing ring is interpreted as the result of a head-on collision with an 
intruder that took place some 1.4 x 10^ years ago. The charac- 
teristics of G„ rule out this galaxy as the intruder. 

a3) A chain of Ho' emitting regions kinematically associated 
with Gj at distances of about 12 kpc have been detected. These 
regions are identified as young, tidally-induced star-forming re- 
gions (some may be candidate tidal dwarf galaxies) formed as a 
consequence of a past interaction. 

B) The VIMOS data of IRAS F121 15-4656 show the po- 
tential of IFS in the kinematic study of the different phases of 
the interstellar gas: the neutral gas traced by the Nal absorption 
doublet, and the ionized gas traced by the Ha or [Nil] emission 
lines. The main conclusions for this system are: 

bl) The neutral and ionized gas morphology show a different 
structure indicating that they have a different spatial distribution. 
The neutral gas closely follows the stellar distribution while the 
ionized gas traces mostly the location of HII regions in the cir- 
cumnuclear and external regions of the galaxy. 

b2) The ionized gas shows typical kinematic characteristics 
of massive rotating disks, with a large amplitude (peak-to-peak 
of 520 km s"'), a maximum of the velocity dispersion at the 
kinematic/light centre, and a major kinematic axis coincident 
with the stellar photometric axis. A dynamical mass of 8.7 x 
10'" Mq is measured inside a region of 6 kpc radius. The neutral 
gas, although consistent on a large scale with a rotating system, 
shows clear departures from pure rotation and it is dynamically 
hotter than the ionized gas, being more affected by random mo- 
tions. 

b3) Despite being dominated by rotation the ionized gas 
shows large-scale (10 kpc), low-amplitude (20-30 km s ') non- 
rotational motion, which spatially correlates with regions of in- 
creased velocity dispersion. The neutral gas also exhibits an even 
clearer spatial correlation between departures from rotation and 
high velocity dispersion values. 
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Table 1. VIMOS Luminous and Ultraluminous Infrared Galaxy sample 
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Double system (dlst.^17") 
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01:22:20.9 
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2 


Kim et al. (2{K)2) /Irregular with prominent tidal tails in WFPC2 image 
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NGC 0633. Double system (dist. ~ 1') 


F02021-2104 
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Kim et al. (2(K)2) / Edge-on spiral in WFPC2 ? 


F02587-6336t 




02:59:43.4 


-63:25:00 


0.264800 


79385 


1345.5 
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Compact Ir in WFPC2 image 


F043 15-0840 


NGC 1614 


04:33:59.8 


-08:34:44 
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69.1 
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11.69 


2 


Arp 186. Spiral with large scale tidal structures 


F05020-2941t 


Abell 3307 


05:04:00.7 
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46168 
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2 


Distorted source with short tidal tails according to Veilleux et al. (2002) 
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Compact Ir in WFPC2 image 


F05 189-2524 
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Tidal structures in the ESO-Danish-1.54m images 
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0.079465 


23823 


360.7 


1501 


12.26 


1 


Double system (dist. ^ 10") in WFC2 image 


F06076-2139 




06:09:45.8 


-21:40:24 


0.037446 


11226 


165.0 


743 


11.67 


1 


On-axis merger 


F06206-6315 




06:21:01.2 


-63:17:23 


0.092441 


27713 


423.3 


1720 


12.27 


2 


Irregular with a prominent arm / tidal tail in WFPC2 image 


F06259-4708 


ESQ 255-IG 007 


06:27:22.0 


-47: 10:54 


0.038790 


11629 


171.1 


769 


11.91 


1 


Triple System 


F06295-1735 


ESQ 557-G002 


06:31:47.2 


-17:37:17 


0.021298 


6385 


92.7 


431 


11.27 





Spiral 


F06592-6313 




06:59:40.2 


-63:17:53 


0.022956 


6882 


100.0 


464 


11.22 





Spiral (Ir), wide pair 


F07027-601 1 


AM 0702-601 


07:03:24.1 


-60:15:23 


0.031322 


9390 


137.4 


626 


11.64 





Spiral 


F07160-6215 


NGC 2369 


07:16:37.7 


-62:20:37 


0.010807 


3240 


46.7 


221 


11.16 





Spiral with irregular inner parts 


F08355-4944 




08:37:01.8 


-49:54:30 


0.025898 


7764 


113.1 


521 


11.60 


2 


Irregular morphology in ACS image 


F08424-3130 


ESQ 432-IG006 


08:44:27.6 


-31:41:50 


0.016165 


4846 


70.1 


329 


11.04 


1 


Double system (otf-set 25") 


F08520-6850 


ESQ 60-IG016 


08:52:29.9 


-69:01:58 


0.046315 


13885 


205.4 


909 


11.83 


1 


Secondary nucleus at — 2". 


F09022-3615 




09:04:12.7 


-36:27:01 


0.059641 


17880 


267.0 


1153 


12.32 


2 


Faint tidal tail in DSS image '? 


F09437+O317 


IC-563/ 564 


09:46:20.6 


+03:03:30 


0.020467 


6136 


89.0 


415 


11.21 


1 


Arp 303. Double system (dist. ~ 90") 


F10015-0614 


NGC-3110 


10:04:02.1 


-06:28:29 


0.016858 


5054 


73.1 


343 


11.31 





Spiral. Interacting with galaxy at SW ? 


F10038-3338 


IC2545 


10:06:35.0 


-33:51:30 


0.034100 


10223 


149.9 


679 


11.77 


2 


Prominent spiral arm or tidal tail in ACS image 


F10173+0828 




10:20:00.2 


+08:13:34 


0.049087 


14716 


218.1 


961 


11.86 


2 


Linear structure in ACS image with evidence of a faint tilted tidal tail 


F10257-4338 


NGC 3256 


10:27:51.3 


-43:54:14 


0.009354 


2804 


40.4 


192 


11.69 


2 


Long tidal tails. Only iimner parts in FoV 


F10409-4556 


ESQ 264-G036 


10:43:07.7 


-46: 12:45 


0.021011 


6299 


91.4 


425 


11.26 





Spiral 


F10485-1447t 




10:51:03.1 


-15:03:22 


0133000 


39872 


625.6 


2363 


12.32 


1 


Kim et al. (2002). Double/triple system 


F10567-4310 


ESQ 264-G057 


10:59:01.8 


- 43:26:26 


0.017199 


5156 


74.6 


350 


11.07 





Spiral 


Fl 1254-4120 


ESO319-G022 


11:27:54.1 


-41:36:52 


0.016351 


4902 


70.9 


333 


11.04 





Spiral. Inner ring. Barred 


F11506-3851 


ESQ 320-G030 


11:53:11.7 


-39:07:49 


0.010781 


3232 


46.6 


221 


11.30 





Spiral 


F12042-3140 


ESQ 440-IG 058 


12:06:51.9 


-31:56:54 


0.023203 


6956 


101.1 


468 


11.37 


1 


Double system. Tidal structures 


F121 15-4656 


ESQ 267-G030 


12:14:12.8 


-47: 13:43 


0.018489 


5543 


80.3 


375 


11.11 





Spiral, wide pair 


F12116-5615 




12:14:22.1 


-56:32:33 


0.027102 


8125 


118.5 


545 


11.61 


2 


Inner spiral-hke structure and outer elliptical morphology in ACS image. 


F12596-1529* 


MCG-02-33-098 


13:02:19.7 


-15:46:03 


0.015921 


4773 


69.0 


324 


11.07 


2 


Sc pec. Close interacting pair 


F 1300 1-2339 


ESQ 507-G070 


13:02:52.3 


-23:55:18 


0.021702 


6506 


94.5 


439 


11.48 


2 


Itregulai; Tidal structures 


F13229-2934 * 


NGC 5135 


13:25:44.0 


-29:50:01 


0.013693 


4105 


59.3 


280 


11.29 





SB spiral. In group 


F14544-4255 -*- 


IC4518 


14:57:42.9 


-43:07:54 


0.015728 


4715 


68.2 


320 


11.11 


1 


Sc pec. close interacting pair 


F17138-1017 * 




17:16:35.8 


-10:20:39 


0.017335 


5197 


75.2 


352 


11.41 


2 


Irregular 


F18093-5744 * 


IC 4687/4686 


18:13:39.6 


-57:43:31 


0.017345 


5200 


75.3 


353 


11.57 


1 


Spiral pec. Close interacting pair 


F21453-3511* 


NGC 7130 


21:48:19.5 


-34:57:05 


0016151 


4842 


70.0 


329 


11.41 


2 


Sa pec ? 


F22132-3705* 


IC5179 


22:16:09.1 


-36:50:37 


0011415 


3422 


49.3 


234 


11.22 





SA(rs)bc, isolated 



" As given in the NASA Extragalactic Database (NED) 

* Luminosity distances using the calculator by Wright (2006) for a lambda-cosmology with Ho=70 kms^'Mpc"', Q.m=03, and £^^=0.7. 
Logarithm of the infrared luminosity, L,>=L(8-1000yum), in units of solar bolometric luminosity, calculated using the fluxes in the four 
IRAS bands as given in Sanders et al. (2003), when available. Otherwise from the IRAS fluxes given in IRAS Point Source and Faint Source 
catalogues (Moshir et al. 1990). 

Morphology class according to a simplified classification scheme similar to that given by Veilleux et al. (2002) for ULIRGs. Note that 
some objects are just single isolated objects without traces of undergoing a merging process, and they are classified as class 0. Galaxy mergers 
in a pre-coalencense phase are classified as class 1, and objects with peculiar morphologies suggesting a post-coalescence merging phase are 
classified as class 2. Note that this classification is based on existing imaging data and in some cases it may change after complete reduction 
and analysis of the sample. 

This system was classified by Veilleux et al. (2002) as a late-stage merger. However, we consider that its morphology can be simply 
interpreted as a nearly edge-on disc galaxy. 

t At the redshift of this object the standard instrumental setting was such that it excluded the Ho- region, but it includes the HyS and the [OIII] lines. 
★ Galaxies with pending observations. 
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Table 2. Main physical properties of IRAS F06076-2139 and IRAS F121 15 -4656 derived from VIMOS IFS 



Quantity 


IRAS F06076-2139 
G„ 


IRAS F06076-2139 
G, 


IRAS F121 15-4656 


V u 


11236 ± 10 


11750 ± 10 


5449 ± 1 1 


^ nucleus,Ha^ 


11186 ± 25 


11737 ± 25 


5453 ± 15 


^ sys,NaD^ 


— 


— 


5460 ± 18 


\T d 


— 


— 


5446 ± 15 




76 ± 5 


61 ± 3 


251 ± 2 


^nuclei! 5^ 


104 ± 10 


63 + 10 


111 ± 10 




5.6 


10.1 


6.0 


rZ 


2.76 


2.73 


7.7 




0.75 


1.7 


8.7 


Mj ' 


5.2 (0.4) 


1.9 (0.15) 


16.5 (1.2) 




1.71 


1.26 


31.3 


T / 


0.56 


0.41 


2.42 


T "I 


4.68 




1.29 




0.4 


0.3 


1.9 


SFR;£ 


81 




22 


log([01]/Hail 
log([NII]/HQ') " 


-1.42 


-1.48 


-1.31 


-0.12 


-0.29 


-0.15 


log([SlI]/Ho') " 


-0.50 


-0.46 


-0.39 



" systemic velocity in km s ' as obtained from the best fit model to the observed Ha velocity field. The errors include both wavelength 
calibration and model fitting uncertainties 

* observed systemic velocity in km s"' as obtained from the nuclear Ho- line 

systemic velocity in km s"' as obtained from the best fit model to the observed NaD velocity field in IRAS F121 15-4656. Uncertainties 
include both wavelength calibration and model fitting 

observed systemic velocity in km s"' as obtained from the nuclear NaD line 
'' rotational velocity in km s"' as obtained from best fit model to the observed Ho- velocity field 
^ velocity dispersion in km s"' measured using the nuclear Ha line profile 

* galactocentric radius in kpc for the rotational velocity given above 

'' efi^ective radius in kpc from the available ACS F814W image for IRAS F06076-2139, and from Lauberts & Valentijn 1989 for IRAS 
F121 15-4656 

' dynamical mass in IO'^Mq as derived from the rotational velocity and galactocentric radius given above 

' dynamical mass in IO'^Mq derived from the effective radius and velocity dispersion values given above. The values in parenthesis are given 
in units of m. (1.4 x 10"Me, Cole et al. 2001) 

* observed flux in 10"'* erg s"' cm"- not corrected for internal extinction 
' in 10"" erg s"' from observed Ha fluxes given above 

in 10" Lo 

" Star formation rate in Mq yr"' assuming SFR = 7.9 x lO"''^ x L^,, 

° Star formation rate in Mq yr"' assuming SFR = 4.5 x 10"'*'' x L;« (Kennicutt 1998). The value for IRAS F06076-2139 represents the 
contribution from the two galaxies, G„ and G, 

not corrected for internal extinction, ratios represent the values obtained for the nucleus 



